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I. INTRODUCTION 



The determination of the classical ground states of interacting many-particle systems 
(minimum energy configurations) is a subject of ongoing investigation in condensed-matter 
physics and materials science While such results are readily produced by slowly 

freezing liquids in experiments and computer simulations, our theoretical understanding of 
classical ground states is far from complete. Much of the progress to rigorously identify 
ground states for given interactions has been for lattice models, primarily in one dimension 
2|. The solutions in d-dimensional Euclidean space M'^ for d > 2 are considerably more 
challenging. For example, the ground state(s) for the well-known Lennard- Jones potential 
in or M.^ are not known rigorously sl. Recently, a "collective-coordinate" approach has 
been employed to study and ascertain ground states in two and three dimensions for a 

n n n 

certain class of interactions [4, |5[. A surprising conclusion of Ref. |^ is that there exist 
nontrivial disordered classical ground states without any long-range order 6|, in addition to 
the expected periodic ones. Despite these advances, new theoretical tools are required to 
make further progress in our understanding of classical ground states. 

In a recent Letter, we derived duality relations for a certain class of soft pair potentials 
that can be applied to classical ground states whether they are disordered or not [3]. Soft 
interactions are considered because, as we will see, they are easier to treat theoretically 
and possess great importance in soft-matter systems, such as colloids, microemulsions, and 
polymers jsH]. These duality relations link the energy of configurations associated with a 
real-space pair potential v{r) to the energy associated with the dual (Fourier-transformed) 
potential. Duality relations are useful because they enable one to use information about the 
ground states of certain soft short-ranged potentials to draw conclusions about the nature 
of the ground states of long-ranged potentials and vice versa. The duality relations also lead 
to bounds on the zero-temperature energies in density intervals of phase coexistence. 

In the present paper, we amplify and extend the results of Ref. j?]. We also study in detail 
a one- dimensional system that Torquato and Stillinger claimed to possess an infinite number 
of structural phase transitions from Bravais to non-Bravais lattices at T = as the density is 
changed . A general set of potential functions that are self-similar under Fourier transform 
are described and studied. We also derive the generalizations of the duality relations for 
three-body as well as higher-order interactions. 
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FIG. 1: (Color online) Left panel: Portion of a Bravais lattice with one particle per fundamental 
cell. Particles are situated on the vertices of the four rhombic fundamental cells shown. Right 
panel: Portion of a periodic non-Bravais lattice with five particles per fundamental cell. 

II. DEFINITIONS AND PRELIMINARIES 

A point process in R'^ is a distribution of an infinite number of points at number density 
p (number of points per unit volume) with configuration ri,r2,...; see Ref. for a 

precise mathematical definition. It is characterized by a countably infinite set of n-particle 
generic probability density functions p„(ri, . . . , r„), which are proportional to the probability 
densities of finding collections of n particles in volume elements near the positions ri, . . . , r„. 
For a general point process, it is convenient to introduce the n-particle correlation functions 
gn, which are defined by 



Since pn = p" for a completely uncorrelated point process, it follows that deviations of 
Qn from unity provide a measure of the correlations between points in a point process. Of 
particular interest is the pair correlation function, which for a translationally invariant point 
process of density p can be written as 



Closely related to the pair correlation function is the total correlation function, denoted by 
h] it is derived from g2 via the equation 




(1) 




(2) 




(3) 
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Since g2{r) — t- 1 as r — t- +oo (r = |r|) for translationally invariant systems without long- 
range order, it follows that h{r) — )■ in this limit, meaning that h is generally an function, 
and its Fourier transform is well-defined. 

It is common in statistical mechanics when passing to reciprocal space to consider the 
associated structure factor S, which for a translationally invariant system is defined by 

Sik) = l + ph{k), (4) 

where h is the Fourier transform of the total correlation function, p is the number density, 
and k = |k| is the magnitude of the reciprocal variable to r. The li- dimensional Fourier 
transform of any integrable radial function /(r) is 

/» = (2.)^/^ fr-VW^^yS^^r, (5) 



Jo " ' {krf'^^-' 

and the inverse transform of f{k) is given by 

Here k is the wavenumber (reciprocal variable) and Jy{x) is the Bessel function of order v. 

A special point process of central interest in this paper is a lattice. A lattice A in M*^ is 
a subgroup consisting of the integer linear combinations of vectors that constitute a basis 



for M , i.e., the lattice vectors p; see Ref. [13| for details. In a lattice A, the space M can 



be geometrically divided into identical regions F called fundamental cells, each of which 
corresponds to just one point as in Figure [H In the physical sciences, a lattice is equivalent 
to a Bravais lattice. Unless otherwise stated, for this situation we will use the term lattice. 
Every lattice has a dual (or reciprocal) lattice A* in which the sites of that lattice are 
specified by the dual (reciprocal) lattice vectors q ■ p = 27rm, where m = ±1, ±2, ±3 ■ • • . 
The dual fundamental cell F* has volume vp* = {2'kY/vf, where vp is the volume of the 
fundamental cell of the original lattice A, implying that the respective densities p and p* of 
the real and dual lattices are related by pp* = l/(27r)'^. A periodic point process, or non- 
Bravais lattice, is a more general notion than a lattice because it is is obtained by placing 
a fixed configuration of points (where A^ > 1) within one fundamental cell of a lattice A, 
which is then periodically replicated (see Figure [I]). Thus, the point process is still periodic 
under translations by A, but the A^ points can occur anywhere in the chosen fundamental 
cell. Although generally a non-Bravais lattice does not have a dual, certain periodic point 



patterns are known to possess formally dual non-Bravais lattices. Roughly speaking, two 
non-Bravais lattices are formal duals of each other if their average pair sums (total energies 
per particle) obey the same relationship as Poisson summation for Bravais lattices for all 
admissible pair interactions; for further details, the reader is referred to [l^ . 



III. DUALITY RELATIONS 



A. Pair Potentials 

For a configuration r'^ = ri, r2, . . . , rjv of ^ 1 particles in a bounded volume y C M*^ 
with stable pairwise interactions, the many-body function 

i=ij=i 

is twice the total potential energy per particle [plus the "self-energy" "^(0)], where v{rij) 
is a radial pair potential function and Vij = \rj — rj|. A pair interaction v{r) is stable 
provided that J2iLi — all > 1 and all G M"^. A nonnegative Fourier 

transform v{k) implies stability, but this is a stronger condition than the former {isij. A 
classical ground-state configuration (structure) within V is one that minimizes f/(r^). Since 
we will allow for disordered ground states, then we consider the general ensemble setting 
that enables us to treat both disordered as well as ordered configurations. The ensemble 
average of U for a statistically homogeneous and isotropic system in the thermodynamic 
limit is given by 

(f/(r^)) =v{r = 0)+p I v{r)g2ir)dr, (8) 

where p = limAr^ooy^oo N/V is the number density and g2{r) is the pair correlation function. 
In what follows, we consider those stable radial pair potentials v{r) that are bounded and 
absolutely integrable. We call such functions admissible pair potentials. Therefore, the 
corresponding Fourier transform v{k) exists, which we also take to be admissible, and 

{U{r^)) = v{r = 0) + pv{k = 0) + p I v{r)h{r)dr. (9) 

JR'i 

Lemma. For any ergodic configuration in R'^, the following duality relation holds: 

/ v{r)h{r)dY = [ v{k)h{k)dk (10) 
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If such a configuration is a ground state, then the left and right sides of (fTOj) are minimized. 
Proof: We assume ergodicity, i.e., the macroscopic properties of any single configuration in 
the thermodynamic hmit A^, — t- +oo with p = N/V = constant are equal to their ensemble- 
average counterparts. The identity flTUl) follows from Plancherel's theorem, assuming that 
h{k) exists. It follows from and fllOl) that both sides of f llOp are minimized for any ground- 
state structure, although the duality relation ( ITOj) apphes to general (i.e., non-ground-state) 
structures. 
Remarks: 

1. The general duality relation f llOl) does not seem to have been noticed or exploited 



before, although it was used for a specific pair interaction in Ref. [13|. The reason 
for this perhaps is due to the fact that one is commonly interested in the total energy 
or, equivalently, the integral of ([8]) for which Plancherel's theorem cannot be applied 
because the Fourier transform of g2{r) does not exist. 



2. 



t is important to recognize that whereas h{r) always characterizes a point process 



12|, its Fourier transform h{k) is generally not the total correlation function of a 
point process in reciprocal space. It is when h{r) characterizes a Bravais lattice A (a 
special point process) that h{k) is the total correlation function of a point process, 
namely the reciprocal Bravais lattice A*. 

3. The ensemble- averaged structure factor is related to the collective density variable 
X]^i6xp(ik ■ Tj) via the expression limjv^.oo(;^| Z^jLi^^Pl^k ■ rj)p) = (27r)'^p5(k) + 
S{k). 

4. On account of the "uncertainty principle" for Fourier pairs, the duality relation ( !T0|) 
provides a computationally fast and efficient way of computing energies per particle of 
configurations for a non-localized (long-ranged) potential, say v{r), by evaluating the 
equivalent integral in reciprocal space for the corresponding localized (compact) dual 
potential v{k). 

Theorem 1. If an admissible pair potential v{r) has a Bravais lattice A ground-state 
structure at number density p, then we have the following duality relation for the minimum 

Umin of U : 

v{r = 0) + Y^v{r) = pv{k = 0) + p ^ 'v{k), (11) 

rGA keA* 



where the prime on the sum denotes that the zero vector should be omitted, A* denotes 



the reciprocal Bravais lattice 16|, and v{k) is the dual pair potential, which automatically 
satisfies the stability condition, and therefore is admissible. Moreover, the minimum Umin 
of U for any ground-state structure of the dual potential v{k), is bounded from above by 
the corresponding real-space minimized quantity Umin or, equivalently, the right side of (ITT]) . 



I.e., 

Umin < Umin = pv{k = 0) + p 2^ v{k). (12) 

keA* 

Whenever the reciprocal lattice A* at reciprocal lattice density p* = p^^(27r)^^ is a ground 
state of v{k), the inequality in ( |T2l) becomes an equality. On the other hand, if an admissible 
dual potential v{k) has a Bravais lattice A* at number density p*, then 

Umin < Umin = P*v{r = 0) + p*Y^v{r), (13) 

reA 

where equality is achieved when the real-space ground state is the lattice A reciprocal to A*. 
Proof: The radially averaged total correlation function for a Bravais lattice, which we now 
assume to be a ground-state structure, is given by 

^"^1'^'^ n=l 

where si(r) is the surface area of a (i-dimensional sphere of radius r, Zn is the coordination 
number (number of points) at the radial distance r„, and 5{r) is a radial Dirac delta function. 
Substitution of this expression and the corresponding one for h{k) into f ITU]) yields 

v{r = Q) + Y,Znv{rn) = pv{k = + pY,Znv{kn), (15) 

71=1 n=l 

where is the coordination number in the reciprocal lattice at the radial distance Rec- 
ognizing that Ylin=i Znv{rn) = Yl'reA^i''") [cQual to twice the minimized energy per particle 
Umin given by ([7]) at its minimum in the limit — )■ oo] a nd Y^ ^-^ Znv{kn) = Yl'keA* "^l^) 



yields the duality relation f lTT]) . The fact that v{r) is stable 15| means that the dual poten- 
tial v{k) is stable since the left side of fll5p is nothing more than the sum given in Ref. 15| 
in the limit — oo, which must be nonnegative. However, the minimum Umin is generally 
not equal to the corresponding minimum Umin associated with the ground state of the dual 
potential v{k), i.e., there may be periodic structures that have lower energy than the recip- 
rocal lattice so that Umin < Umin- To prove this point, notice that U for any non-Bravais 



lattice by definition obeys the inequality Umin < U. However, because the corresponding 
Fourier transform h{k) of total correlation function h{r) of the non-Bravais lattice in real 
space generally does not correspond to a point process in reciprocal space (see Remark 2 
under Lemma 1), we cannot eliminate the possibilities that there are non-Bravais lattices in 
reciprocal space with If lower than Umin- Therefore, the inequality of (11 2p holds in general 
with equality applying whenever the ground state structure for the dual potential v{k) is the 
Bravais lattice A* at density p*. Inequality (|T3l) follows in the same manner as f|T2|) when 
the ground state of the dual potential is known to be a Bravais lattice. 
Remarks: 

1. Whenever equality in relation ( !T2|) is achieved, then a ground state structure of the 
dual potential v{k = r) evaluated at the real-space variable r is the Bravais lattice A* 
at density p* = p"^(27r)^'^. 

2. The zero- vector contributions on both sides of the duality relation (|TT|) are crucial in 
order to establish a relationship between the real- and reciprocal-space "lattice" sums 
indicated therein. To emphasize this point, consider in the well-known Yukawa 
(screened-Coloumb) potential v{r) = exp{—K,r)/r, which has the dual potential v{k) = 
(47r)/(K^ -|- k"^). At first glance, this potential would seem to be allowable because the 
real-space lattice sum, given on the left side of f lTTj) . is convergent. However, the 
reciprocal-space lattice sum on the right side does not converge. This nonconvergence 
arises because v{r = 0) is unbounded. Equality of "infinities" is established, but of 
course this is of no practical value and is the reason why we demand that an admissible 
potential be bounded. 

3. Can one identify specific circumstances in which the strict inequalities in (I12p and (II 3p 
apply? In addition to the theorem below that provides one such affirmative answer to 
this question, we will also subsequently give a specific one- dimensional example with 
unusual properties. 

Theorem 2. Suppose that for admissible potentials there exists a range of densities over 
which the ground states are side by side coexistence of two distinct structures whose parent- 
age are two different Bravais lattices, then the strict inequalities in f lT2|) and f|T3|) apply at 
any density in this density-coexistence interval. 
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Proof: This follows immediately from the Maxwell double-tangent construction in the U- 
plane, which ensures that the energy per particle in the coexistence region at density p 
is lower than either of the two Bravais lattices. 

As we will see, the duality relations of Theorem 1 will enable one to use information 
about ground states of short-ranged potentials to draw new conclusions about the nature of 
the ground states of long-ranged potentials and vice versa. Moreover, inequalities (fT2|) and 
f|T3|) provide a computational tool to estimate ground-state energies or eliminate candidate 
ground-state structures as obtained by annealing in Monte Carlo and molecular dynamics 
simulations. In the ensuing discussion, we will examine the ground states of several classes 
of admissible functions, focusing under what conditions the equalities or strict inequalities 
of the duality relations f lT^ and ffT^ apply. 

B. Three-Body and Higher-Order Interactions 

The aforementioned analysis can be extended to establish duality relations for many- 
particle systems interacting via three-body and higher-order interactions. For simplicity of 
exposition, we begin with a detailed construction of the three-body duahty relations and 
then generalize to the higher-order case. 

We consider a statistically homogeneous A^-particle interaction $jv(r^) with one-, two-, 
and three-body contributions fi, V2{rij), and v^irij^ra^rji), respectively. With the conven- 
tion that V2{rij) = Vsi^Vij, rij, 0) and Vi = V2{0) = t'3(0, 0, 0), we may write 

*iv(r^) = ^ X1^3(rij,ri^,rj7), (16) 

where ^3 is symmetric, bounded, and short-ranged. Taking the ensemble average of this 
function implies 

($Ar(r^)) = ui 3p / g2{v)v2{r)dv + j g:i{v,s)v^{r, s,\t - s\)drds, (17) 

involving averages over single particles, pairs, and triads. Duality relations for the former 
two contributions have already been considered, and we therefore direct our attention to the 
last term in (IT7|) . 

Since g3{r,s) — > 1 as |r|, |s|, and |r — s| — )■ 00, this function is generally not integrable, 
and we therefore introduce the associated three-body total correlation function /i3(r,s) = 

9 



(73 (r, s) — 1. Application of a double Fourier transform and Plancherel's theorem implies the 
following three-body analog of the Lemma ffTOj) : 

/ h3{r,s)v3{r,s,\r-s\)drds= . . \. ■ / h3{k,q)v3{k,q)dkdq, (18) 

where 

/(k,q)=/ e.p(->k.r-.q.s)/(r,sMris. (19) 

One can verify directly that the following relationship defines the three-body correlation 
function for any statistically homogeneous A^-particle point pattern: 

p'gsir, s) = / i 5^ 5(r - r,,)5(s - r,,) \ . (20) 

For a Bravais lattice, ergodicity should hold, and we can re-write (1201) as 

p2(73(r,s) = ^'5(r-r,)5(s-r,), (21) 

where the set {fj} in the summations includes all points of the lattice excluding the origin. 

The dual Bravais lattice will possess a three-particle correlation function of the form 
^3(k, q) = 1 + p^/;,3(k, q), where p is the real space number density. Substituting fl2T|) and 
the corresponding for the dual Bravais lattice into (JTSj) gives the following duality relation 
for three-particle interactions: 

1^1+^ '"■i{rj,r^, \\Yj -r^ll) = p^vi + {;3(k^,k„), (22) 

where we have defined vi = ^3(0, 0). 

The extension of this analysis to higher-order interactions is straightforward. Specifically, 
we consider a n-particle bounded, symmetric, and short-ranged potential fn('"i2, • • • , ''"in) 
with a statistically homogeneous point distribution and the associated Plancherel identity 

/i„(ri2, . . . , ri„)f„(ri2, . . . , ri„)rfri2 ■ ■ ■ drin = 

'n-l)d 

{ — ] / /in(ki,...,k„_i)t;„(ki,...,k„_i)(iki---(ik„_i. 

(23) 

The n-particle correlation function of a Bravais lattice is 

•••,rn_i)= Yl '^(ri-raJ---5(r„_i-r„„_J, (24) 



{a}n-l 
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where denotes all sets of n — 1 distinct vectors in a Bravais lattice, excluding the 

origin, and a indexes the lattice points. Using this relationship, we find the following general 
n-particle duality relation: 



{"In-l 



^^1 + ^ t;„({ka}) 



(25) 



where v\ = f„({0}). 



IV. APPLICATIONS 

A. Admissible functions with compact support 

Recently, the ground states have been studied corresponding to a certain class of oscillat- 
ing real-space potentials v{r) as defined by the family of Fourier transforms with compact 
support such that v{k) is positive for < k < K and zero otherwise . Clearly, v{k) is an 
admissible pair potential. Siito [5] showed that in three dimensions the corresponding real- 
space potential f (r), which oscillates about zero, has the body-centered cubic (bcc) lattice 
as its unique ground state at the real-space density p = 1/{8^/2tc^) (where we have taken 
K = 1). Moreover, he demonstrated that for densities greater than 1/(8\/2-it^), the ground 
states are degenerate such that the face-centered cubic (fee), simple hexagonal (sh), and 
simple cubic (sc) lattices are ground states at and above the respective densities l/(6-\/37r^), 
v^/(16V27r3), and l/(8V27r3). 

The long-range behavior of the real-space oscillating potential v{r) might be regarded to 
be unrealistic by some. However, since all of the aforementioned ground states are Bravais 
lattices, the duality relation ( ITTl) can be applied here to infer the ground states of real-space 
potentials with compact support. Specifically, application of the duality theorem in and 
Siito's results enables us to conclude that for the real-space potential v{r) that is positive 
for < r < D and zero otherwise, the fee lattice (dual of the bcc lattice) is the unique 
ground state at the density a/2 and the ground states are degenerate such that the bcc, sh 
and sc lattices are ground states at and below the respective densities (3a/3)/4, 2/y/3, and 
1 (taking D = 1). Specific examples of such real-space potentials, for which the ground 
states are not rigorously known, include the "square- mound" potential [v{r) = e > Q 
for < r < 1 and zero otherwise] and what we call here the "overlap" potential, which 

11 



corresponds to the intersection volume of two d-dimensional spheres of diameter D whose 
centers are separated by a distance r, divided by the volume of a sphere. The latter potential, 
which has support in the interval [0, D), remarkably arises in the consideration of the variance 
in the number of points within a spherical "window" of diameter D for point patterns 
in and its minimizer is an open problem in number theory [l^. The (i-dimensional 
Fourier transforms of the square mound and overlap potentials are e2'^^'^ Jd/2{k) / (kTcY^"^ and 
2'^7r'^/2r(l + d/2)J^^^{k/2)/k'^, respectively, with D = 1. Figure |2] shows the real-space 
and dual potentials for these examples in three dimensions. The densities at which the 
aforementioned lattices are ground state structures are easily understood by appealing to 
either the square-mound or overlap potential. The fee lattice is the unique ground state 
at the density a/2 because at this value (where the nearest-neighbor distance is unity) and 
lower densities the lattice energy is zero. At a slightly higher density, each of the 12 nearest 
neighbors contributes an amount of e to the lattice energy. At densities lower than \/2, 
there is an uncountably infinite number of degenerate ground states. This includes the 
bcc, sh and sc lattices, which join in as minimum-energy configurations at and below the 
respective densities (3v^)/4, 2/y/3, and 1 because those are the threshold values at which 
these structures have lattice energies that change discontinuously from some positive value 
(determined by nearest neighbors only) to zero. Moreover, any structure, periodic or not, 
in which the nearest-neighbor distance is greater than unity is a ground state. 

However, at densities corresponding to nearest-neighbor distances that are less than unity, 
rigorous prediction of the possible ground-state structures is considerably more difficult. For 
example, it has been argued in Ref. 10[ (with good reason) that real-space potentials whose 
Fourier transforms oscillate about zero will exhibit polymorphic crystal phases in which the 
particles that comprise a cluster sit on top of each other. The square-mound potential is 
a special case of this class of potentials and the fact that it is a simple piecewise constant 
function allows for a rigorous analysis of the clustered ground states for densities in which 
the nearest-neighbor distances are less than the distance at which the discontinuity in v{r) 
occurs [10 1 . 
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FIG. 2: (Color online) Left panel: The square-mound potential v{r) = G(l — r) and the three- 
dimensional overlap potential v{r) = 0(1 — r) [l — 3r/2 -|- , where Q{x) is the Heaviside step 
function. Right panel: Corresponding dual potentials v{k) = vr'^/^ J3/2(^)/(2A;)^^^ (square-mound; 
scaled by vr^/G for clarity) and v{k) = Gvr^ [j3/2(A;/2)] ^ //c^ (overlap). 

B. Nonnegative admissible functions 

Another interesting class of admissible functions are those in which both v{r) and v{k) 
are nonnegative (i.e., purely repulsive) for their entire domains. The "overlap" potential 
discussed above is an example. 



1. One- Dimensional Overlap Potential 

Here we examine the one-dimensional ground-state structures associated with the dual 
of the so-called overlap potential 

vir)=(l-^)Q{D-r), (26) 

which is equal to the intersection volume, scaled by D, of two rods of radius D/2 with 
centers separated by a distance r. The dual potential is 

'sm{kD/2y " 



v{k) = D 



(27) 



{kD/2) 

Figure |3] shows that both potentials are bounded and repulsive. However, while the overlap 
potential possesses the compact support [0,/^], the dual potential is long-ranged with a 
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FIG. 3: Left panel: One-dimensional overlap potential v{r) = 0(1 — r)(l — r). Right panel: 
Corresponding dual potential v{k) = 4sin^(A;/2)//c^. 



countably infinite number of global minima determined by the zeros k* = 2m7r/Z} (m G N) 
of sin{kD/2). Torquato and Stillinger have shown 13|] that the unique ground state of the 
d = 1 overlap potential is the integer lattice with density p = l/D; Theorem 1 therefore 
implies that the integer lattice at reciprocal density p* = D/{2-k) is the unique ground 
state of the dual potential ( 1271) . This result intuitively corresponds to placing each point 
in an energy minimum of the dual potential, thereby driving the total potential energy 
to zero. This argument immediately implies that the integer lattice at reciprocal density 
p* = D / (27rm) for all m G N is also a ground state of the dual potential; however, the 
ground states at intermediate densities are generally non-Bravais lattices and have heretofore 
been unexplored. Based on these observations, previous work has suggested that the dual 
interaction ( |27j) undergoes an infinite number of structural phase transitions from Bravais 
or simple non-Bravais lattices to complex non-Bravais lattices over the entire density range 

Q. 

We have characterized the ground states of the dual overlap potential numerically using 



the MINOP algorithm 



19j, which applies a dogleg strategy using a gradient direction when 



one is far from the energy minimum, a quasi- Newton direction when one is close, and a linear 
combination of the two when one is at intermediate distances from a solution. The MINOP 
algorithm has been shown to provide more reliable results than gradient-based algorithms for 
similar many-body energy minimization problems 20|. We fix the length L of the simulation 
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FIG. 4: (Color online) Illustrative portions of numerically-determined ground state configurations 
of the dual potential v{k) (p7|) with D = 2tt and densities p = I (upper left), p = 1/2 (upper right), 
p = 2/3 (lower left), and p = 4/3 (lower right). The particles have been given a small but finite 
size for visual clarity. Note that the p = 4/3 configuration is a "clustered" integer lattice with 
more than one particle occupying certain lattice sites. 

box and use a modified version of the dual potential 



where is the number of particles. Note that L/A provides the unit of length for the 
problem, allowing us to control the density of the resulting configuration by varying A. 

For the case A = 1, we have numerically verified that the the integer lattice is the unique 
ground state (up to translation) of the dual potential; indeed, direct calculation shows that 
the integer lattice minimizes the potential energy (1281) for all A G N as expected from the 
arguments above. However, we have also identified degenerate ground states that are non- 
Bravais lattices; these systems are shown in Figure |H Our results suggest that for A > 1 the 
ground states are complex superpositions of Bravais lattices with a minimum inter-particle 
spacing determined by A. Additionally, the conjectured infinite structural phase transitions 
are not observed in this density range owing to the high degeneracy of the ground state. We 
remark that although the integer lattice is a ground state for any A G N, it is never observed 
in our numerical simulations because the energy landscape possesses a large number of global 
minima. Furthermore, although the ground states for integral and non-integral values of A 
are visually similar, we emphasize that the integer lattice is never a ground-state candidate 
for A ^ N. For A < 1, the ground states are more difficult to resolve numerically because 
finite-size effects become more pronounced in this region; justification for this behavior is 
provided below. Nevertheless, we observe a "clustered" integer lattice structure in which 
several points occupy a single lattice site for 0.5 < A < 1. 

Our numerical results suggest an exact approach to characterizing the ground states of 




(28) 
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the dual potential (127|) . For simplicity and without loss of generality, we will henceforth 
consider the scaled pair interaction 

Hr)-{^)\ (29) 

corresponding to a normalized dual potential with D = To facilitate the approach to the 
thermodynamic limit, we first examine a compact subset of M subject to periodic boundary 
conditions. The entropy of this system for p < 1 can be determined by relating the problem 
to the classic model of distributing balls into M > N jars such that no more than one 
ball occupies each jar (Fermi-Dirac statistics). Specifically, choosing the parameter A > 1 in 
(128|) is equivalent to choosing a density p = l/A<lin the general problem fl29|) . Therefore, 
for any A > 1, there are M = AN "jars" for the particles ("balls"). Assuming that 
the particles are indistinguishable, the number of distinct ways of distributing the particles 
into the M potential energy minima is the binomial coefficient (^) = (^). For A^ large 
(approaching the thermodynamic limit), Stirling's formula implies that the entropy S is 

= AATln^^^^^ +A^ln(A-l), (31) 

where we have chosen units with k-^ = 1. Rearranging terms and substituting A = 1/p for 
the density, we find 

S/N = p-' \n{p-') - ip-' - 1) ln(p-i - 1), (32) 

which is fixed in the thermodynamic limit and is plotted in Figure |5l Note that S/N \ as 
p 1, which is expected from the observation that the integer lattice is the unique ground 
state at unit density. This unusual residual entropy refiects the increasing degeneracy of the 
ground state with decreasing density and implies that in general the aforementioned infinite 
structural phase transitions from Bravais to non-Bravais lattices are not thermodynamically 
observed. Instead, one finds an increasing number of countable coexisting ground-state 
structures as seen in our numerical energy minimizations. 

For p > 1, determination of the ground states of the dual potential ( 129|) is nontrivial since 
it is no longer possible to distribute all of the particles into potential energy wells. Therefore, 
Fermi-Dirac statistics are no longer applicable for the many-particle system. Nevertheless, 
we can make some quantitative observations concerning the ground states in this density 

16 




FIG. 5: (Color online) Entropy per particle S/N as a function of density p for the dual potential 
^(A;) in ^ with D = 2-n. 



regime. We first consider the scenario of adding one particle to a local region, subject to 
periodic boundary conditions, of the integer lattice of unit spacing. Since the potential 
energy minima of the pair interaction ( l29l) occur on the sites of the integer lattice, the 
total potential energy cannot be driven to its global minimum. Symmetry of the lattice 
implies that, without loss of generality, we can limit the location ^ of the particle to the 
interval [0,0.5]. Since the energy of the underlying integer lattice is zero and the particle, 
by construction, will not interact with periodic images of itself, the total potential energy 
of the system after addition of the particle is exactly 



^ ( sin[7r(C + n)] y ^ [ sin[7r(l - j + n)] 



n=0 

sin^l TT. 



<i + n) ; ■ ^ V 7r(l-e + n) 



(33) 
(34) 



where ^l)'^^]ix) is the trigamma function 2l|. From the reflection property of the trigamma 
function |2l|], the latter expression is exactly equal to unity [= "5(0)] for any value of the 
parameter ^. 

Determination of the ground state then depends on "relaxing" the system by making a 
small perturbation 7 < 1 in the underlying integer lattice (see Figure [6]). The energy E' of 
this perturbed system is then parametrized by the displacements 7 and ^ as in Figure [6] and 
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FIG. 6: (Color online) Upper: Schematic illustrating the "relaxation" of the integer lattice toward 
the ground state of the dual potential (j29p for p slightly greater than unity. Lower: Energy 
landscape associated with local perturbations of the integer lattice. 

is given by 

^, ^ / sin[7r(l - 7 + n)] \ f sin[7r(l + 7 + n)] \ f sin[7r(l -j + n) ] ^ ^ 

7r(l-7 + n) ) f-V vr(l + 7 + n) ) ^\ 7^(1-^ + ^) 



n=0 ^ ^ ' ' n=0 ^ ^ ' ' n=0 

^ / sin[7r(l + 7 + n)] y ^ / sin[7r(7 + ^ + n)] V 

vr(l + 7 + n) y V 7r(7 + ^ + n) / ^ 



n=0 



2 + 5(7 + -^^(7)-^5(0, (36) 



where we have utilized reflection and recurrence relations for the polygamma function [21 1 
with v{r) given by ( l29l) . Figure [6] illustrates that E' possesses a unique minimum value 
= 0.777216 at (7,^ = (0.279376,0.279376). Because the pair interaction (El) is long- 
ranged, it is unclear if -E^in can be further decreased by additional local deformation of the 
structure. Nevertheless, our analysis suggests that the ground state structures at densities 
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slightly above unity are perturbed integer lattices with "defects" in the crystal structure. 

Upon reaching p = 2, it is possible to "stack" two integer lattices with unit spacing for an 
energy per particle E/N = 1/2; interestingly, the long-range nature of the pair potential fl29p 
implies that these lattices can be mechanically decoupled from each other without increasing 
the energy of the system. Furthermore, the symmetry of this configuration implies that no 
local perturbation of the lattice structure can decrease the energy per particle, meaning that 
this "stacked" integer lattice and its translates within [0, 1] are at least local minima of the 
pair interaction ( 129|) : a similar argument will hold for any p G N. The energy of this stacked 
configuration is 



E/N={a-l){a-2)A/2 + {a-l)[l-{a-l)A] A = 1/p = 1/a (a G N). (37) 



1. If the stacked integer lattices are global minima of the pair interaction ( l29ll for any 
number density p G N, then the ground states are unique (up to translation of layers) at 
these densities, and the residual entropy will therefore vanish. However, for p ^ N there 
is a combinatorial degeneracy associated with local deformations of the underlying 
integer lattice, implying that the residual entropy is nonanalytic over the full density 
range. This behavior in combination with the thermodynamic relation 



where /3 is the thermal expansion coefficient and kt is the isothermal compressibility, 
suggests that there exist densities where the ground state exhibits negative thermal 
expansion as T — )■ 0. 

2. One special case of the aforementioned "stacked" integer lattice configurations occurs 
when multiple particles occupy the same lattice sites (i.e., with no translation between 
layers). For these "clustered" integer lattices, pair interactions are localized to include 
only those particles on the same lattice site, meaning that there are no long-range 
interactions for these systems. However, we have seen that the inclusion of long-range 
pair interactions, such as with the p = 2 integer lattice, does not affect the total energy 
of the system. Since relative displacements between layers are uniformly distributed 
on [0,1], the average displacement of 0.5 indeed corresponds to the p = 2 integer 
lattice. 



Remarks: 




(38) 
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3. Our results imply that numerical methods are in general not appropriate for identifying 
the ground states for p > 1 since truncation of the summation fl33l) (e.g., with the 
minimal image convention) breaks the translational degeneracy of the system. 

4. The ground states of the overlap potential v{r) = (1 — r)0(l — r) also exhibit rich 
behavior for p > 1. Since the interactions are localized to nearest neighbors, one can 
verify that addition of a particle to the unit density integer lattice increases the energy 
of the system by one unit, regardless of the position of the particle. However, unlike 
the dual potential fl2^ . no local perturbation of the integer lattice can drive the system 
to lower energy, resulting in a large number of degenerate structures. 

5. The two-dimensional ground states of the generalized dual overlap potential 

d{r) = A7T[Ji{kD/2)/kf (39) 



have also been numerically investigated 



22| : the topology of the plane significantly 



increases the difficulty in analytically characterizing the ground-state configurations. 



2. Gaussian-Core Potential 

Another interesting example of nonnegative admissible functions is the Gaussian core 

n 

potential v{r) = eexp[— (r/cr)^] [23], which has been used to model interactions in polymers 



9|, |2J|. The corresponding dual potentials are self-similar Gaussian functions for any d. The 
potential function pairs for the case d = 3 with e = 1 and a = 1 are v{r) = exp(— r^) 



and v{k) = tt"'/^ exp(— A;^/4). It is known from simulations 23| that at sufficiently low 
densities in M.^, the fee lattices are the ground state structures for v{r). It is also known 



that for the range < p < tt^^^^, fee is favored over bcc 25]. If equality in f|T2l) is achieved 
for this density range, the duality theorem would imply that the bcc lattices in the range 
(47r)^'^/^ < p* < oo (i.e., high densities) are the ground state structures for the dual potential. 
Lattice-sum calculations and the aforementioned simulations for the Gaussian core potential 
have verified that this is indeed the case, except in a narrow density interval of fcc-bcc 
coexistence 0.17941 < p < 0.17977 around p = vr"'^/^ ^ 0.17959. In the coexistence interval, 
however, the corollary states the strict inequalities in (1121) and (fT3|) must apply. Importantly, 
the ground states here are not only non-Bravais lattices, they are not even periodic. The 
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ground states are side-by- side coexistence of two macroscopic regions, but their shapes and 
relative orientations are expected to be rather comphcated functions of density, because they 
depend on the surface eriergies of grain boundaries between the contacting crystal domains. 
Proposition 9.6 of Ref. [263 enables us to conclude that the integer lattices are the ground 
states of the Gaussian core potential for all densities in one dimension. Note that in R^, 
the triangular lattices apparently are the ground states for the Gaussian core potential at 
all densities (even if there is no proof of such a conclusion), and therefore would not exhibit 
a phase transition. Similar behavior has also been observed in four and eight dimensions, 
wWe .elf-dua, E, lattice, a. the apparent g.ou„d .ate. Q,Q. Co., K..a.. 

and Schiirmann have recently identified non-Bravais lattices in five and seven dimensions 
with lower ground-state energies than the densest known Bravais lattices and their duals 
in these dimensions IJ]. Interestingly, these non-Bravais lattices, which are deformations 
of the D'^ and packings, possess the unusual property of formal self- duality, meaning 
that their average pair sums (total energies per particle) obey the same relation as Poisson 
summation for Bravais lattices for all admissible pair interactions. It is indeed an open 
problem to explain why formally-dual ground states exist for this pair potential. 

It is also instructive to apply our higher-order duality relations f l25|) to the simple example 
of a three-body generalization of the aforementioned Gaussian-core potential. Specifically, 
we consider a three-body potential of the form 

^3(^12, ^^23) = exp(-r^2 - - r^g) = exp [-2(r^2 + r^g - na • rig)] . (40) 

Applying a double Fourier transform to this function shows that the dual potential, given 
by 

^3(k, q) = j exp [-{e + + l|k - q|n/12] , (41) 

is self-similar to fHOj) . As with the two-body version of the Gaussian-core potential, this self- 
similarity implies that if a Bravais lattice is the ground state of the three-body Gaussian-core 
interaction at low density, then its dual lattice will be the ground state at high density with 
the exception of a narrow interval of coexistence around the self-dual density p* = (S/vr^)'^/^. 
However, we have been unable to find either numerical or analytical studies of the ground 
states of this higher-order interaction in the literature, and determining whether it shares 
ground states with its two-body counterpart is an open problem. 
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C. Completely monotonic admissible functions 



A radial function /(r) is completely monotonic if it possesses derivatives for all 

n > and if (— l)"/*^")(r) > 0. A radial function /(r) is completely monotonic if and 
only if it is the Laplace transform of a finite nonnegative Borel measure /i on [0, oo], i.e., 
/(^) ~ /o°° ^~^*'^/^(^) 28|. Not all completely monotonic functions are admissible (e.g., the 
pure power-law potential in M"^ is inadmissible). Examples of completely monotonic 
admissible functions in include exp(— ar) for a > and l/(r + a)^ for a > 0, P > d. 
Importantly, the Fourier transform f{k) of a completely monotonic radial function /(r) is 
completely monotonic in k"^ 29|. 

Remarkably, the ground states of the pure exponential potential have not been inves- 
tigated. Here we apply the duality relations to the real-space potential v{r) = exp(— r) 
in M'^ and its corresponding dual potential v{k) = c{d)/(l + A;2)('^+i)/2 [where c{d) = 
2''"ii^^~^^^'^T{{d + l)/2)], which has a slow power-law decay of l/k'^^^ for large k. Note 
that the dual potential is a completely monotonic admissible function in k"^, and both 
v{r) and v{k) also fall within the class of nonnegative admissible functions. We have per- 
formed lattice-sum calculations for the exponential potential for a variety of Bravais and 
non-Bravais lattices in and M^. In R^, we found that the triangular lattices are fa- 
vored at all densities (as is true for the Gaussian core potential). If equality in f fT2|) is 
achieved, then the triangular lattices are also the ground states for the slowly decaying dual 
potential v{k) = 2tc/{1 + k"^)^^"^ at all densities. In R'^, we found that the fee lattices are 
favored at low densities (0 < p < 0.017470) and bcc lattices are favored at high densities 
(0.017470 < p < oo). The Maxwell double-tangent construction reveals that there is a 
very narrow density interval 0.017469 < p < 0.017471 of fcc-bcc coexistence. We see that 
qualitatively the exponential potential appears to behave like the Gaussian core potential. 
If equality in f lT2|) applies outside the coexistence interval, then the duality theorem would 
predict that the ground states of the slowly-decaying dual potential v{k) = Svr/ (1 + A;^)^ are 
the fee lattices for < p* < 0.230750 and the bcc lattice for 0.230777 < p* < oo. Note 



that in one dimension, it also follows from the work of Cohn and Kumar [26| that since the 
integer lattices are the ground states of the Gaussian potential, then these unique Bravais 
lattices are the ground states of both the exponential potential and its dual evaluated at 
k = r (i.e., v{r) = 2/(1 + r^)). 
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Cohn and Kumar 



26| have rigorously proved that certain configurations of points inter- 



acting with completely monotonic potentials on the surface of the unit sphere in arbitrary 
dimension were energy-minimizing. They also studied ways to possibly generalize their 
results for compact spaces to Euclidean spaces and conjectured that the densest Bravais 
lattices in R'^ for the special cases d = 2, 8 and 24 are the unique energy-minimizing con- 
figurations for completely monotonic functions. These particular lattices are self-dual and 
therefore phase transitions between different lattices is not possible. Note that if the ground 



states for completely monotonic functions of squared distance in M (the Gaussian 



unction 



26| that 



being a special case) can be proved for any d > 2, it immediately follows from Ref. 
the completely monotonic functions of distance share the same ground states. Thus, proofs 
for the Gaussian core potential automatically apply to the exponential potential as well as 
its dual (i.e., f (r) = c{d)/{l + r^)('^+^)/^) because the latter is also completely monotonic in 

^2 



Based upon the work of Cohn and Kumar [26|, it was conjectured that the Gaussian 
core potential, exponential potential, the dual of the exponential potential, and any other 
admissible potential function that is completely monotonic in distance or squared distance 
share the same ground-state structures in for 2 < d < 8 and d = 24, albeit not at the 
same densities [TJ. Moreover, it was also conjectured for any such potential function, the 
ground states are the Bravais lattices corresponding to the densest known sphere packings 
30| for < p < pi and the corresponding reciprocal Bravais lattices for p2 < P < cxd, 
where pi and p2 are the density limits of phase coexistence of the low- and high-density 
phases, respectively. In instances in which the Bravais and reciprocal lattices are self-dual 
{d = 2,4,8 and 24) pi = p2, otherwise p2 > pi (which occurs for d = 3,5,6 and 7). The 
second conjecture was recently shown by Cohn and Kumar to be violated for d = 5 and 
d = 7. 



V. SELF-DUAL FAMILIES OF PAIR POTENTIALS 

Our discussion of the Gaussian core model above suggests that one can exactly map the 
energy of a lattice at density p to that of its dual lattice at reciprocal density p* for pair 
potentials that are self-similar (defined below) under Fourier transform. Here we provide 
additional examples of self-similar pair potentials, including radial functions that are eigen- 
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functions of the Fourier transform. Only some of these results are known in the mathematics 
literature 3l| , and this material has not previously been examined in the context of duality 
relations for classical ground states. 



A. Eigenfunctions of the Fourier transform 

Pair potentials that are eigenfunctions of the Fourier transform are unique in the con- 
text of the duality relations above since they preserve length scales for all densities; i.e., 
v{k) = Xv{k) with no scaling factor /i in the argument. We therefore briefly review these 
eigenfunctions and the associated eigenvalues for radial Fourier transforms. In order to 
simply the discussion, we will adopt a unitary convention for the Fourier transform in this 
section 

/(k) = J exp(-zk ■ r)/(r)^;r = ^{/}(k), (42) 

which differs from our previous usage only by a scaling factor. The slight change in notation 
(/ instead of /) is intended to clarify which convention is being used. 

The eigenfunctions of the Fourier transform for d = 1 can be derived from the generating 
function for the Hermite polynomials, which, when scaled by a Gaussian, is given by 

exp (—X /2 + 2tx -t^) = ^i — ] exp (— Hn{x). (43) 

Taking the Fourier transform of both sides, one obtains 

(^i-^ ^ exp (-t^) £exp {-(l/2)[x2 - x(4t - 2ik)]} dx = J2{j^)^ {exp(-xV2)//„(x)} 

(44) 

implying 

exp(-A;V2) [^^j Hnik) = ('^'j ^ {exp(-xV2)i7„(x)} . (45) 

By collecting powers of t in ( 145|) . we immediately conclude 

^?{exp(-a;V2)i7„(x)} = (-^^ exp(-A;V2)i/„(A;), (46) 

thereby identifying both the eigenfunctions and eigenvalues of the d = 1 Fourier transform. 
Note that the eigenvalues are real when n is even. 
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We now seek eigenfunctions of the radially-symmetric Fourier transform, defined here as 



f{r)r 



d-l 



J, 



(kr) 



dr 



(47) 



for an isotropic function /(r). Direct substitution shows that /(r) = exp(— r^/2) is an 
eigenfunction for all d with eigenvalue 1. Other eigenfunctions of the Fourier transform can 
be identified by noting that they are also eigenfunctions of the c?- dimensional Schrodinger 
equation for the radial harmonic oscillator 



, , . fd-l\ d , . . 



dr 



dr 



where we have used the relation 



fd-1 

+ 



dr 



d 
dr 



(48) 



(49) 



for radially-isotropic functions in d dimensions. The eigenvalues of the Schrodinger equation 
are En = n + d/2 for some n G N U {0}. The general solutions to ( HHi) are then given by 

Ci{d) exp (— r^/2) L^^^'^~^\r'^) {k = n/2 for n even) 



(50) 



where (x) is the associated Laguerre polynomial [2l| and Ci (d) is a dimension- dependent 
constant. Note that for = 1 



(j)k{x) = exp (-r^/2) L 



(-1/2), 



r oc 



exp (-rV2) H2ki 



(51) 



and we recover the even d = 1 eigenfunctions of the harmonic oscillator. 

To determine the eigenvalues of the radial Fourier transform, we note that if / is an 
eigenfunction, then it must be true that 



f{k) = cf{k) 

for some eigenvalue c. However, it is also true that 

f(k) = kk)=ci\k) = ^f{k). 



(52) 



(53) 



Equation fl53|) implies that either c = ±1 or f{k) = 0; for a nontrivial solution we conclude 
that the eigenvalues of the radially-symmetric Fourier transform are ±1, which is in contrast 
to the general case on Mf^. This result is exactly consistent with the constraint that the index 
n of an eigenstate of the radial Schrodinger equation ( HHl) be even. Note that when c = — 1, 
the Fourier transform changes the nature of the interaction (i.e., repulsive to attractive and 
vice- versa) . 
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B. Poly-Gaussian potential 



The results above can be extended to include linear combinations of eigenfunctions of the 
Fourier transform; furthermore, we can generalize these functions to be simply self-similar 
under Fourier transform, meaning that length scales are not preserved by the transformation. 
Specifically, our interest is in functions for which: 

v{k) = Xviixk), (54) 

where A and fi are constants. 

As an example, we consider the Gaussian pair potential of the Gaussian core model 

f{r,a)=exp{-{r/af). (55) 

The corresponding Fourier transform is given by 

f{k,a) = / exp(— (r/cr)^) exp(ik ■ r) 

= n'/^a'^expi-a^e/A). (56) 

Now consider a pair potential v{r) that is a linear combination of two Gaussians as follows: 

v{r)=A,f{r,a,) + A2f{r,a2). (57) 

Its Fourier transform is 

d{k) = 7T^/^ [Aiaf exp{-aie/A) + A2al exp{-alk^ / A)] . (58) 

In order for v{r) to be self-similar under Fourier transformation, the constants /z and A 
of f l5^ must satisfy the following two equations for all x: 

a\exp{~a\x^ / A) = AA2 exp(— (/ix/(72)^) (59) 
/A) = \Aiexp{-{fix/aif). (60) 

These equations will be satisfied by requiring 



= A = (2.^P, A,--0j,, (61) 



leaving three independent parameters: /i, cxi, and Ai. 
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The example extends to any even number of Gaussian components. Let 



2n 



v{r) = J^exp(-(r/(Tj)2), 



(62) 



where the a,- are ordered by magnitude: 



< (Ji < (J2 < • • • < o"2n < +00 



(63) 



The corresponding Fourier transform is given by 



2n 



i){k) = vr^/^ J2 Ajo] exp(-ajA;V4). 



(64) 



In order to ensure self-similarity, the terms can be paired and subject to the relations of 
the type fl6Tl) . On account of the ordering condition fl63|) . we pair terms with indices j and 
2n — j + 1, 1 < j < n, and hence require 



C"2r!,-j+l 



2/i 

a,- 



(65) 



It is also possible to include an additional Gaussian to make an odd number in total. 
This additional term must effectively pair with itself, so that 



(To = (2/i) 



1/2 



(66) 



where the corresponding parameter Aq is uncontrained and the subscript refers to the 
"odd"' term. 

These relations suggest an extension to the case of a continuous distribution of Gaussian 
widths as follows: 



v{r) = / A{a) 
Jo 



^2 ^3 
exp( -) + ^^^^^expi 



a^' (2/i)3/2 
The corresponding Fourier transform is given by 

,21,2 ^3 



da. 



(67) 



v{k) = {27Tfif^ / A{a) 
Jo 



exp( —) + -——exp{ 



a2 ^ (2/x)3/2 



4/x2 



da 



(68) 



as required for self-similarity, where A = (27r/i) 



3/2 
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VI. DISCUSSION AND CONCLUSIONS 



In this work we have derived duahty relations for interactions of arbitrarily high order 
that can be applied to help quantify and identify classical ground states for admissible po- 
tentials that arise in soft-matter systems. We have applied the duality relations for different 
classes of admissible potential functions, including potentials with compact support, nonneg- 
ative functions, and completely monotonic potentials. Among these classes, the completely 
monotonic functions offer a new category of potentials for which the ground states might 
be identified rigorously. In particular, we seek a proof of the conjecture that functions in 
this class share the same ground-state structures in R'' for 2 < d < 8 and d = 24, albeit 
not at the same densities. No counterexample for this conjecture has been found to date. It 
should also be emphasized that the examples of admissible functions examined here are by 
no means complete. 

We have also identified a set of pair potentials on the line related to the overlap function 
that exhibit a "stacking" phenomenon at certain densities in the ground state. This behavior 
leads to an unusual mechanical decoupling between layers of integer lattices due entirely to 
the form of the interaction. These systems, previously thought to exhibit an infinite number 
of structural phase transitions from Bravais to non-Bravais structures [t], likely possess rich 
thermodynamic properties such as negative thermal expansion as T — )■ 0. Since overlap 
potentials ..se in a variety of contexts, inc.ud.ng the covering and quantize, problems Q 
and the identification and design of hyperuniform point patterns [13|], further studies of these 
systems are warranted. 

Toward this end, we plan to explore whether analogous duality relations can be estab- 
lished for positive but small temperatures by studying the properties of the phonon spectra 
of admissible potentials. The development of such relations would provide a unique and 
useful guide for mapping the phase diagrams of many-particle interactions, including those 
functions belonging to the class of "self-similar" potentials that we have introduced here. 
Indeed, with the exception of the Gaussian core model 271], little is known about the ground 
states and phase behaviors of self-similar functions. Since most of these potentials contain 
both repulsive and attractive components, these interactions have direct implications for 
spatially inhomogeneous solvent compositions that simultaneously induce repulsion and at- 
traction among macromolecules in solution. We expect that as the methodology continues 
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to develop, duality relations of the type we have discussed here will play an invaluable role 
in understanding these complex physical systems. 
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